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We confront the Two-Higgs-Doublet Model with a variety of experimental constraints as well as 
theoretical consistency conditions. The most constraining data are the B Xs'y decay rate (at low 
values of Af^±), and Ap (at both low and high Mu±). We also take into account the BB oscillation 
rate and i?;,, or the width T{Z bb) (both of which restrict the model at low values of tan /3), and 
the B~ tVt decay rate, which restricts the model at high tan/3 and low M[j±. Furthermore, the 
LEP2 non-discovery of a light, neutral Higgs boson is considered, as well as the muon anomalous 
magnetic moment. Since perturbative unitarity excludes high values of tan (5, the model turns out to 
be very constrained. We outline the remaining allowed regions in the tan /3-Af^± plane for different 
values of the masses of the two lightest neutral Higgs bosons, and describe some of their properties. 



I. INTRODUCTION AND NOTATION 

As compared with the Standard Model (SM), the Two- 
Higgs-Doublet Model (2HDM) allows for an additional 
mechanism for CP violation ^1,, ^, ^, 4]. This is welcome, 
in view of baryogenesis Q , and one of the main reasons 
for continued interest in the model. 

Several experimental constraints restrict it. The B — B 
oscillations and branching ratio Ri, exclude low values of 
tan /?, whereas the B — > rate excludes low values 
of the charged-Higgs mass, Mh± ■ The precise measure- 
ments at LEP of the p parameter constrain the mass 
sphtting in the Higgs sector, and force the masses to be 
not far from the Z mass scale Q. These individual con- 
straints are all well-known, but we are not aware of any 
dedicated attempt to combine them, other than those of 
. The present study aims to go beyond that of 9] , by 
using more complete and more up-to-date experimental 
results, as well as more accurate theoretical predictions 
for the above quantities. 

From the theoretical point of view, there are also vari- 
ous consistency conditions. The potential has to be posi- 
tive for large values of the fields El . We also require 
the tree-level Higgs-Higgs scattering amplitudes to be 
unitary [l^, [H, Il4| . Together, these constraints dramat- 
ically reduce the allowed parameter space of the model. 

The present study is limited to the 2HDM (H), which 
is defined by having one Higgs doublet ($2) couple to 
the up- type quarks, and the other ($1) to the down-type 
quarks [15| . 

We write the general 2HDM potential as: 



y = ^($t$,)2 + ^($t$,)2 



A3($I$i)($t$2) 
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A4($I$2)($^*l) 

i{m?i(ci>l<f>i) + 



A5($I$2)' + h.c. 



,($l$2)+ll.C. 



-fm22 



(1.1) 

(1>5*2)} 



Thus, the Z2 symmetry will be respected by the quar- 
tic terms, and Flavour-Changing Neutral Currents are 
constrained [ll|. We shall refer to this model (without 
the Ae and A7 terms) as the 2HDM5. The more general 
model, with also Ae and A7 couplings, will be discussed 
elsewhere. 

We allow for CP violation, i.e., A5 and 771^2 may be 
complex. Thus, the neutral sector will be governed by a 
3x3 mixing matrix, parametrized in terms of the angles 
Qfi, a2 and as in (4. [l7|: 



Cl C2 Sl C2 S2 

i? = I -(Ci S2 S3 + S1 C3) C1C3-S1S2S3 C2S3 

-Cl S2 C3-I-S1 S3 -(Ci S3-I-S1 S2 C3) C2C3^ 

where ci = cosai, si = sinai, etc., and 



(1.2) 



(In ref. [H, the angles are denoted as a — ai, a}, — 
ct2, etc = 0(3.) For a discussion of this parameter space, 
including the CP-nonviolating limits, see [Ti|. We will 
use the terminology "general 2HDM" as a reminder that 
CP violation is allowed. The present study extends that 
of 9] also in this respect. 

Rather than taking the parameters of the potential 
(jl.ip to describe the model, we take the two lightest neu- 
tral Higgs boson masses, Mi and M2, together with the 
mixing angles (|1.3p . the charged Higgs boson mass, M^± , 
and tan/? as our basic parameters. (The third neutral 
Higgs boson mass, M3, is then a derived quantity.) 

It is convenient to split the constraints into three cat- 
egories: 

(i) Theoretical consistency constraints: positivity of 
the potential [Tol . [Tlj and perturbative unitarity 
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(ii) Experimental constraints on the charged-Higgs sec- 
tor. These all come from B-physics, and are due to 
b sj, B-B oscillations, and B rvr- They are 
all independent of the neutral sector. 

(iii) Experimental constraints on the neutral sector. 
These are predominantly due to the precise mea- 
surements of non-observation of a neutral Higgs 
boson at LEP2, Ap, and = j{g ~ 2)^. 

The first and third categories of constraints will depend 
on the neutral sector, i.e., the neutral Higgs masses and 
the mixing matrix. The second category is due to phys- 
ical effects of the charged-Higgs Yukawa coupling in the 
B-physics sector. These are "general" in the sense that 
they do not depend on the spectrum of neutral Higgs 
bosons, i.e., they do not depend on the mixing (and pos- 
sible CP violation) in the neutral sector. 

When considering the different experimental con- 
straints, our basic approach will be that they are all 
in agreement with the Standard Model, and simply let 
the experimental or theoretical uncertainty restrict pos- 
sible 2HDM contributions (this procedure yields lower 
bounds on the charged-Higgs mass, possibly also other 
constraints). An alternative approach would be to actu- 
ally fit the 2HDM to the data. This will not be discussed 
in the present paper. 

This paper is organized as follows. In Sec. [IT] we dis- 
cuss the general constraints of positivity of the Higgs po- 
tential, together with tree-level unitarity of Higgs-Higgs 
scattering amplitudes. The impact of these constraints 
is displayed in the tan/3-M^± plane for a few repre- 
sentative values of neutral Higgs boson masses. Next, 
in Sec. Illli we discuss the constraints coming from the 
-B-physics experiments, in particular the B Xs"/, 
B Tv-T and B-B oscillations. Sec. IIVI is devoted to 
various experimental constraints that depend on details 
of the neutral-Higgs sector. In Sec.|Vl we combine all the 
constraints and in Sec. IVII we give some characteristics 
of the surviving parameter space. In Sec. IVIII we specu- 
late on possible future experimental constraints, and then 
summarize in Sec. I Villi 

II. GENERAL THEORY CONSTRAINTS 



loss of detailed information can be compensated for by 
the "overview" obtained by the following procedure: 

1. Pick a set of neutral- Higgs-boson masses, (Mi, M^) 
together with /i^. 

2. Scan an = ni x n2 x 713 grid in the ai-a2-a^ 
space, and count the number j of these points that 
give a viable model. (Alternatively, one could scan 
over A'^ random points in this space.) 

3. The ratio 

Q = j/N, < g < 1, (2.2) 

is then a figure of merit, a measure of "how al- 
lowed" the point is, in the tan/3-7\f^± plane. If 
Q = 0, no sampled point in the a — (ai,Q!2,a3) 
space is allowed, if Q = 1, they are all allowed. An 
alternative measure 

Q+=j/N+, Q+>Q, (2.3) 

counts in the denominator only those points A^+ for 
which positivity is satisfied. 

Of course the 2HDM, if realized in nature, would only 
exist at one point in this parameter space. However, we 
think the above quantities Q and (3+ give meaningful 
measures of how "likely" different parameters are. 



A. Reference masses 

We shall impose the conditions of positivity, unitar- 
ity and experimental constraints on the model, for the 
different "reference" mass sets given in Table U (and vari- 
ations around these). For each of these mass sets we scan 
the model properties in the a = (ai,Q;2,a3) space. (All 
scans have been performed over a 200 x 200 x 100 grid. 
We note that other scanning procedures might be more 
efficient [l8|.) From these reference masses, some trends 
will emerge, allowing us to draw more general conclu- 
sions. 



In the general CP-non-conserving case, the neutral sec- 
tor is conveniently described by the three mixing angles, 
together with two masses (Mi, M2), tan/3 and 



Remi2 
2 cos /3 sin /3 ' 



(2.1) 



We shall here project these constraints from the 
multi-dimensional parameter space onto the tan P-Mjj± 
plane. Such a projection of information from a multi- 
dimensional space onto a point in the tan (3-Mjj± plane 
can be done in a variety of ways, all of which will lead 
to some loss of information. However, we feel that this 



Name 


Mi[GeV] 


MafGeV] 




"100-300" 
"150-300" 
"100-500" 
"150-500" 


100 
150 
100 
150 


300 
300 
500 
500 


[±(200)"] 
[±(200f ] 
[±(200f ] 
[±(200f ] 



TABLE I: Reference masses. 



In the 2HDM5, with ImAs 7^ 0, the two input masses 
will together with a and tan /3 determine M3 [l^j . Spec- 
ifying also Mh± and /x^, all the A's can be determined. 
(For explicit formulas, see (T9|.) 
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B. Positivity and unitarity 

Let us first discuss the effect of imposing positivity. 
Actually, we will use the term "positivity" to refer to 
the non-trivial conditions M| > and M2 < M3 to- 
gether with l/($i,$2) > as |<i>i|, |$2| —> 00. We shall 
henceforth refer to the set of points in the a space where 
positivity is satisfied, as a+ . In Table II of 19] we report 
percentages Q of points in a space for which positivity is 
satisfied. For the mass parameters of Table I, the frac- 
tion is around 30%. For small and negative values of /i^, 
"most" of the exclusion provided by the positivity con- 
straint is due to the conditions M| > and M3 > M2, 
without the explicit conditions on the A's discussed, for 
example, in Appendix A of 11]. 



(M,,lvl;) = (100,500) GeV 



(M,,MJ = (1 50,300) GeV 




1 10 

(M,,M;) = (100,500) GeV 
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tan/? 



FIG. 1: Percentage of points Q+ in the a+ space that satisfy 
unitarity. Four sets of (Mi,M2) values are considered, as 
indicated. The contours show Q+ = 0, 20%, 40% and 60% 
(upper panels) and 0, 5%, 10% and 15% (lower panels). All 
panels: ^ = 0. Yellow region: Q+ > 0. Also shown, are the 
90% and 95% C.L. exclusion contours from Fig. (2] 

In Fig. [Tl for /j? = 0, we study the effects of imposing 
unitarity llilil [13. This has a rather dramatic effect 
at "large" values of tan (3 and Mh± ■ While the general 
constraints on the charged-Higgs sector, to be discussed 
in Sec. IIIIl exclude low values of tan (3 and Mh± , the 
constraints of unitarity exclude high values of these same 
parameters. Only some region in the middle remains not 
excluded. For (Mi,M2) = (100,300) GeV and /x^ = 0, 
unitarity excludes everything above tan /3 ~ 5 (for any 
value of M[j±), and above M^± ~ 650 GeV (for any 
value of tan /?) . 

We shall refer to the set of a values for which unitarity 
as well as positivity are satisfied as d G a+. For M2 — 
300 GeV (upper panels), the percentage of points in a+ 



space for which unitarity is satisfied, reaches (at low tan f3 
and low M[j±) beyond 60%, whereas for M2 — 500 GeV 
(lower panels), it only reaches values of the order of 15- 
20%. 

The domains in which (3+ > depend on /i^: For 
negative values of /x^, the region typically shrinks to lower 
values of tan/3, for positive values of n"^ it extends to 
larger values of tan/3. When /i^ ~ 0, and tan/3 < 5, 
significant fractions of the d space are allowed. However, 
for large positive values of /i^, when also large values 
of tan /3 are allowed, only small domains in d remain 
allowed, as discussed in [1^. They can actually be hard 
to find in a regular, equidistant scan over d. 

The unitarity constraints are conveniently formulated 
in terms of the different weak isospin and hypercharge 
channels [l^ . At large values of tan /3 it turns out to be 
the isospin-zero, hypercharge-zero channel that is most 
constraining. 



III. GENERAL CONSTRAINTS FROM THE H± 
SECTOR 

The precise B-physics experiments provide severe con- 
straints on the charged-Higgs sector, excluding low values 
of Mh± and tan/3. We shall here discuss the three most 
severe constraints of this kind that are independent of 
the neutral sector. 



A. Xs-y 

The B — > Xsj branching rate was early found to con- 
strain the allowed charged-Higgs-boson masses, but also 
to be very sensitive to QCD effects. At leading logarith- 
mic order (LO), it is given by [13, [Ml ill : 



\Vt*sVtb\^ 6ae 



X B{B 



TTg{z) 



(3.1) 



where the first factor is a ratio of CKM matrix ele- 
ments, g{z = m'^/ml) is a phase space factor, and 

c!f'^'^^ (lib) is an effective Wilson coefficient, evaluated 
at the B-meson scale, fit- This effective Wilson coef- 
ficient is obtained from the relevant ones at the elec- 
troweak scale, cf'''(^o), i — 2,7,8 _23;|, where the ef- 
fects of the 2HDM enter. Certain linear combinations 
are denoted "effective" coefficients, they are defined such 
that e.g. C7°''°*(/i) includes all one-loop contributions to 
b S7, i.e., also those of four-quark operators 23j. The 
effective Wilson coefficient was early found to be quite 
sensitive to the scale relevant to i3-meson decay, chang- 
ing by ±25% if the scale ^f, is varied by a factor of two 
in either direction around nib — 5 GeV [H, . At the 
NLO, this scale sensitivity is however significantly re- 
duced M, M M M MmMMM- 
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The additional contributions due to the 2HDM can at 
the weak scale be described by diagrams involving 
exchange, and depend on this mass, Mh± , as well as on 
the Yukawa couplings, i.e., on tan/3. At leading logarith- 
mic order in QCD, they are discussed in [2ll.[34l|. Some of 
these additional terms can be enhanced by factors cot^ /3 
from the Yukawa coupling squared. They all vanish 
linearly in m^/Mfj±, i.e., for Mfj± m^. 

The NLO results for the 2HDM have been studied by 
many authors, see [13, Isil, [sH]. As compared with 
the LO calculation, it has been found that the NLO ef- 
fects weaken the constraints on the allowed region in the 
tan j3-Mu± plane [H, [H, , the bound on Mfj± is sig- 
nificantly relaxed. 

Apart from minor effects, this calculation has for the 
SM been carried to the next-to-next-to-leading order 
(NNLO) [131. At the NNLO, the branching ratio can 
be written as fSU Sal : 



X B{B XceVe)cKp, 



{P{Eo) + NiEo)} 

(3.2) 



where P and N denote perturbative and non-per- 
turbative effects that both depend on the photon lower 
cut-off energy Eq. In the LO limit, the term P reduces to 
the square of the effective Wilson coefficient c!f''°^{fib) 



in (|3.ip . whereas NLO and NNLO contributions include 
effects due to gluon exchange and emission, and require 
the summation over a bilinear expression involving also 
other Wilson coefficients. The factor C accounts for 
mc-dependence associated with the semileptonic decay 
B X^eDe 

The SM prediction of Misiak et al. is (3.15 ± 0.23) x 
10"'' for Eo = 1.6 GeV [13, H^, if aU errors are added 
in quadrature. This is to be compared with the recent 
experimental results, which are averaged to 3.55 x 10^^ 
[4(I|, with an uncertainty of 7-7.5%, again with statis- 
tical and systematic errors added in quadrature. An- 
dersen and Gardi ,41.i| advocate a different approach to 
the resummation of the perturbation series, by "Dressed 
Gluon Exponentiation" , which accounts for multiple and 
soft collinear radiation. These effects are particularly im- 
portant for large photon energies. Their approach yields 
(3.47 ± 0.48) x 10"'', i.e., they find a 14% uncertainty. 
Becher and Neubert also introduce further corrections to 
the calculation of the photon spectrum, and find a rather 
low value, (2.98 ± 0.26) x 10"'* [42], leaving more room 
for new physics. 

The 2HDM contribution is positive, a finite value for 
the charged Higgs mass would thus bring the results of 
[stI . H^I in closer agreement with the experiment. We 
shall however take the attitude that these numbers are 
compatible and compare the uncertainty in the exper- 
imental result and the SM prediction with the 2HDM 
contribution. 

In the NNLO, the perturbative contribution P is ob- 
tained via the following three steps [s^: 



1 . Evaluation of the Wilson coefficients at the "high" 
(electroweak) scale, /^o These coefficients 
are expanded to second order in as and rotated 
to "effective" Wilson coefficients Cf{fio) 

The 2IIDM effects enter at this stage, at lowest 
order in the Wilson coefficients C7(/io) and Cs(po)- 

2. Evaluation of the "running" and mixing of these op- 
erators, from the high scale to the "low" (B-meson) 
scale. This is where the main QCD effects enter 
via a matrix U that is given in terms of powers of 
ry = a,(Mo)/a.(Mh) [H, S H . 



3. Evaluation of matrix elements at the low scale 
[H, 1131, which amounts to constructing P{Eo) of 
Eq. ^ from the Cf{fit)- 

We adopt the scale parameters of (37| : 

Aio = 160 GeV, = 2.5 GeV, = 1-5 GeV. (3.3) 

Actually, our treatment of the higher-order effects has 
been simplified compared to that described in [13, HI] , in 
the sense that: (i) We determine the contribution p^^^^''" 
of [si] by their Eq. (4.10), using results of [H [13, , but 



read off Pr 



(2)rcm 



5 (valid at the "default" scales) from 



their Fig. 2. (ii) Rather than explicitly including the 
0{Vub) and electroweak corrections, we adopt the corre- 
sponding numerical values of -1-1% [U [111 and —3.7% 
[Slj, [3^, respectively. 

Furthermore, while the "matching" in the SM is per- 
formed to second order in Og, we include, in addition to 
the dominant, lowest order, 2HDM effects, also the first- 
order (in asip-o)) contributions. In fact, we include the 
latter at the level of "effective" Wilson coefficients, fol- 
lowing [3^, and take their "matching scale" fiw as /iQ. 
Characteristic relative magnitudes of the LO and NLO 
2HDM contributions are given in Table [Til for tan (3=1 
and 10, and Mh± = 300 GeV and 600 GeV. The LO 
2HDM contribution is measured with respect to the full 
SM value, whereas the NLO 2IIDM contribution (accord- 
ing to |30|) is measured relative to the SM value plus the 
LO 2HDM contribution. The missing 0{a1{^Q)) correc- 
tions appear unimportant. 





tan/3 = 1 


tan/3 = 10 


[GeV] 


LO 


NLO 


LO 


NLO 


600 


18.2 


-3.7 


16.3 


-3.5 


300 


41.8 


-5.5 


36.5 


-5.1 



TABLE II: Relative (in per cent) 2HDM contributions to the 
branching ratio at LO and NLO (Borzumati & Greub). 



As input parameters for (|3.2|1 . in addition to (|3.3|1 . 
we take the CKM ratio to be 0.9676 [ii], B{B 
XccDe) = 0.1061 [53 (see also mt{fj,o) = 162 GeV 
(corresponding to a pole mass of rrit — 171.4 GeV), 
TOe(mc) = 1.224 GeV [H and mf = 4.68 GeV [Hf, 
yielding mc{fJ.c) = 1.131 GeV and mc{^c)l'mb = 0.242. 
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For the non-perturbative part, N{Eo), we follow refs. [38j 
and [HjIIJI- Our result for the branching ratio, Eq. (I3.2p . 
in the SM limit (M^i oo), is 3.12 x lO"^. We col- 
lect in Table Hill the various values obtained by different 
authors. 



Authors 


B 


a 


reference 


HFAG 


3.55 


0.24 ± ■ ■ ■ 


[40] 


M.M. et al. 


3.15 


0.23 


[37] 


J. A. & E. G. 


3.47 


0.48 


m 


T. B. & M. N. 


2.98 


0.26 


[42] 


present work 


3.12 







TABLE III: Branching ratios B{B Xs"f) and uncertainties 
in units of 10"*, all for Ej > Eo ^ 1.6 GeV. 



We define a measure of the amount by which the 
2HDM would violate the agreement with the SM value. 



xl 



(3.4) 



Here, B{B -'^s7)2HDM denotes the 2HDM predic- 
tion. As noted above, it will depend on tan /3 and Mjj± , 
whereas B{B Xs7)rcf denotes a reference value, taken 
to be the averaged experimental value, 3.55 x lO"** 

For the uncertainty that enters in (|3.4p . we adopt the 
value 



a[B{B Xsj)] = 0.35 x 10" 



(3.5) 



which corresponds to the experimental and (SM) theo- 
retical uncertainties [st'I added in quadrature. For the 
2HDM, the theoretical studies have not been carried 
to the same level of precision, but we recall that the 
2HDM-specific NLO-level contributions only modify the 
over-all branching ratios by 0(5%). These numbers give 
xLn = 1-52. 

With the above choice, we can determine exclusion re- 
gions in the tan/3-Mjy± plane. Such regions are shown 
in Fig. [5] for = 4.61 and 5.99, corresponding to confi- 
dence levels of 90% and 95%, respectively, for 2 degrees 
of freedom (tan/3 and Mh±). The obtained bound dif- 
fers from that given by Misiak et al. [IJ even though 
the branching ratio is reproduced to a reasonable preci- 
sion. The reason is as follows. For compatibility with 
the treatment of the other experimental constraints, we 
consider a two-dimensional constraint, whereas they con- 
sider a single-sided one-dimensional constraint.^ 

By adopting the approximate SM-value p^^^-*"^*"" = 5, 
we only include part of the 2HDM-specific contributions 
to P^^\ A more well-defined procedure would be to leave 



600 



> 

O 



400 



200 



(off plot) ^ 




tan/? 



FIG. 2: Excluded regions at small Mff± and small tan/3 due 
to constraints from B — > Xsy and B-B oscillations. Dashed: 
90% C.L., solid: 95% C.L. The corresponding curves for B 
TV are off the plot, to the lower right. Colored: joint exclusion 
at 90% and 95% C.L. 



out all 2HDM-specific contributions to P^^^ ? However, 
that procedure leads to an exclusion limit about 30 GeV 
higher in A/^±, i.e., more of the parameter space would 
be excluded. We have chosen the more conservative ap- 
proach, which leads to less exclusion. 



B. B ^ XrVr and B rVr 

Charged Higgs bosons would contribute to the decay 
h — > ctVt (3.6) 

from which the bound 



tan/3 < 0.52 GeV x Af^± 



(3.7) 



has been obtained [5al at the 2cr level. This bound, which 
would exclude a corner to the lower right outside the 
region shown in Fig. [2] is actually irrelevant, since such 
high values of tan /3 are excluded by the other constraints 
to be discussed in Sec. IIVI 

Charged Higgs-boson exchange also contributes to 



B- 



(3.8) 



which similarly provides a bound of the kind (j3.7p . The 
measurement gives B{B~ tDt = (1.79 ± 0.71) x 10~^ 
[56| . where we have added in quadrature symmetrized 



1 We shall refer to = 4.61 and 5.99 as 90% and 95% C.L., 
respectively, even though this conversion to a probability only ^ We are grateful to M. Misiak for discussions on this and related 
applies to the simple, ideal, case. issues. 
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statistical and systematic errors. With a Standard-Model 
prediction of (1.59 ± 0.40) x lO^'', 



B{B- 



= 1.13 ±0.53. 



(3.9) 



Interpreted in the framework of the 2IIDM, one finds [5 

^2 



rn 2HDM 



1 - 



(3.10) 



We take ms — 5.28 GeV, and formulate a measure as 



\b — yTU 



[(jirH cxp)]^ 



(3.11) 



It follows that two sectors at large values of tan f3 and low 
values of M^± will be excluded. This bound is stronger 
than that of p.7p . It has some relevance in the lower 
right of Fig. [51 when all effects are added (see Sec lIIIDj) . 



B—B oscillations 



The precisely measured Bd-Bd oscillations are at low- 
est order given by the formula (571 



with the Inami-Lim functions 158 



S- 



2HDM 



S 



WW 



25. 



WH 



S 



HH- 



(3.12) 



(3.13) 



The contributions proportional to diagrams with the 
exchange of one or two charged-Higgs-bosons are de- 
noted SwH and Shh- They are proportional to 
cot^/3 and cot^ /?, and vanish when M^± ^ rrit, as 



)\og{milMjj^) 



and ml/Mfj^, respectively. 



At the NLO, the corresponding result has also been ob- 
tained [5^: Bb and rj (denoted 772 in ^5^) receive correc- 
tions of order a^. In particular, the 0{as) contribution 
to rj2 becomes a non-trivial function of tan (3 and Mb± ■ 
A few comments are here in order: (i) The factor of 2 
for the second term in p.l3p has been adopted to follow 
the notation of [5^ . (ii) There is a book-keeping problem 
with the expression for L^'-^) in Eq. (A. 20) of [53|- Since 
the last term in that expression, proportional to a quan- 
tity denoted HH^ has an explicit coefficient l/tan**/?, 
the Shh in Eqs. (A.21) and (A.24) for HH^^ should be 
replaced by tan^ f3 x Shh- (iii) There is a discrepancy in 
a quantity denoted 2WW^l\ between Eq. (A. 16) in ^ 
and the later PhD thesis of the same author. We have 
chosen to take the formula given in the thesis. At the 
level of 772, it amounts to a difference of the order of 2%."^ 



^ We are grateful to A. Buras, U. Jentschura and F. Krauss for 
correspondence on these points (ii) and (iii). 



The Oipts) corrections to 77 introduce a variation of rj 
(or 772) from 0.334 at tan/3 = 0.5 and Mh± = 200 GeV 
to 0.552 at tan/3 — 50. (We have adopted the over- 
all normalization of 77 such that it agrees with the SM 
value 0.552 60] at tan/3 = 50.) However, the product 
rj X S2HDM varies by a factor of 2.7 over this same range, 
as compared with a factor of 4.5 for S2HDM itself. Thus, 
the inclusion of the 0{as) QCD corrections reduce the 
sensitivity of Aran to charged- Higgs contributions. In 
other words, the inclusion of NLO corrections weakens 
the constraints on the 2HDM at low values of tan (3. 

Recently, also the Bg-Bg oscillation parameter AniB^ 
has been measured |61i]. It is given by an expression 
similar to (|3.12p . except that the CKM matrix elements 
are different: 



ArriB. 



Gl 



I ?\yth\^fBBBmBrjMl^ S2HDM , 



67r2' 

(3.14) 

and Bb and 771 b , now referring to Bg , are numerically 
different. 

^ /2 

The quantities Jb BJ that appear in (|3.12|) and (|3.14|) 
are determined from lattice QCD studies, and rather un- 
certain. We shall use the values adopted recently by Ball 
and Fleischer 16(1, 



fB.B],/^ = (0.244 ±0.026) GeV, 
fg^ Bg^^ = (0.295 ± 0.036) GeV, 



(3.15) 



based on unquenched calculations. In fact, these uncer- 
tainties are the dominant ones. 

For each of these observables, one can form a mea- 
sure of the deviation from the SM: 



^B,-B, 



(A771_B^2HDM - AmB^Su)'^ 

WiAmB.W 



(3.16) 



We adopt the attitude that the measurements of A777,b^ 
and AruB^ furnish determinations of the GKM matrix el- 
ements \Vtd\ and |Vts| that are compatible with the SM, 
and simply require that the additional 2HDM contribu- 
tions do not spoil this consistency. With the assumption 
that the uncertainties in (|3.15[) are the dominant ones, 
we have 



cr(A77isJ = 21% X Attib^, 
a{AmBj = 24% x Auib,- 



(3.17) 



However, with the theory error (j3.15p being the domi- 
nant one, we can not claim that the measurements of 
Amsa ^'^'i Attib^ furnish independent constraints on the 
model. Therefore, we consider only the contribution 
from the more constraining one of these two measure- 
ments, namely AmB^ with a 21% uncertainty: 



Xl-B 



-By 



(3.18) 



Contours of X^-s shown in Fig. O As mentioned 
above, they are more generous than the corresponding 
bounds based on the LO theory. 
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D. Combining general constraints from the H 
sector 

In order to quantify the extent to which a particular 
point in the tan (3-AIif± plane is forbidden, we form a 
as follows: 

Xgcneral ~ Xb~*s^ Xb^TV ^ Xb — B' (3.19) 

The subscript "general" refers to the fact that constraints 
depending on the neutral sector are not yet taken into 
account. Yet, no choice of neutral-sector parameters can 
avoid these constraints. 

Figure [2] shows excluded regions in the tan (3-Mh± 
plane, due to the constraints of -B ^ ^s7, B~ —> tD-t and 
B-B oscillations at 90% (dashed) and 95% C.L. (sohd), 
as well as combined (indicated in yellow and red at 90% 
and 95% C.L., respectively). Earlier versions of such ex- 
clusion plots can be found in [1, [3l| . 

IV. EXPERIMENTAL CONSTRAINTS ON THE 
NEUTRAL SECTOR 

We now turn to the constraints coming from experi- 
ments related to the neutral Higgs sector. These are not 
"general" , they will depend on the parameters of this sec- 
tor, namely neutral Higgs masses and mixing angles. We 
will study, as representative cases, those given in TablelH 

A problem with point 2 of the procedure of Sec. |TT1 is 
the following. The conditions of positivity and unitarity 
are absolute'*, with a certain fraction of the points in 
the a+ space satisfying these, whereas the experimental 
constraints are statistical, i.e., a certain parameter point 
may violate an experimental observation by Icr, 2a or 
more. There is no obvious way to combine this informa- 
tion. 

As a first step, we will consider one of these experimen- 
tal constraints on the neutral sector at a time, and display 
the way it excludes parts of the tan (3~Mb± plane by the 
following simple consideration. Let x| be the contribu- 
tion to due to a particular experimental observable 

2 _ (Ci.2HDM — C'i.rcf)^ x 

~ HOiW ' ^ ■ ' 

where the "reference" value Oi^-^ci will be either the ex- 
perimental value or the SM value. 

We shall consider the following observables: i?f,, the 
branching ratio for Z ^bb^ (in Sec lTyT)) : the LEP2 
non-discovery of a light neutral Higgs boson [b^, [S^I 
(in Sec. IIVB[) : p, the LEP determination of the rela- 
tion between the Z and W masses [1^ (in Sec. IIV CP : 



* One might consider models for which perturbativc unitarity is 
not satisfied, but we shall not do that here. 



and = ^{g — 2)^, the precise Brookhaven determina- 
tion of the muon anomalous magnetic moment [66i] (in 
Sec. IIVD[) . Then, the criterion we will consider in this 
first step, can be expressed in terms of the following quan- 
tity: 

• For fixed tan (3 and Mb± take 

X? = min X- (4.2) 

where is minimized over the part a of the a+ 
space for which positivity and also unitarity are 
satisfied. 

The minimization over a finds the point where the model 
can most easily accommodate the particular experimen- 
tal constraint, with the chosen masses (Mi,M2) and 
fi^ fixed. This point in a will in general differ from 
one experimental constraint to another. We will estab- 
lish 90 and 95% C.L. allowed regions in the ia.n [3-Mfj± 
plane, corresponding to one of these observables at a 
time. These will subsequently be combined. However, 
the over-all allowed regions will in general be less than 
the intersection of the individual ones, since the latter 
will correspond to different points in d. 

This procedure preserves none of the probabilistic in- 
formation of Fig. [1] which shows that some region of the 
tan/3~M//± plane contains more possible solutions than 
some other region. All focus is here on the one "best" 
point in a, where x^ is lowest. 

A. Ri, constraint (LEP) 

The 2HDM-specific contributions to Rb are of two 
kinds. At low values of tan/3, the exchange of charged 
Higgs bosons is important, whereas at high values of 
tan/3 the exchange of neutral Higgs bosons is important 
[g^]- For the general CP non-conserving case, the latter 
contributions are given in [ll[. 

For this observable, the reference value in (|4.ip is not 
well defined, since Rb [i^ is part of the electroweak ob- 
servables from which a SM Higgs mass is fitted. Hence, 
we take 

2 ( Ai?f,,2HDM\^ (A Q^ 

where Ai?b^2HDM refers to the 2HDM-specific contribu- 
tions to this quantity [ll|, [63| and with (j{Rb) ~ 0.05% 
the experimental uncertainty [63 |. 

We show in Fig.[3l for two sets of (Mi, M2) values, and 
fi = 0, how this constraint removes a sliver of low-tan /3 
values. The allowed regions are here cut off at tan/3 > 7 
due to the unitarity constraint discussed in Sec. lH Bl The 
neutral-Higgs-exchange contribution to Rb is in this case 
of no importance. 

Next, we show in Fig. [H the corresponding allowed re- 
gions for /i = 200 GeV. The large-tan /3 region is then 
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/i=0 '^b 
(M,,M;) = (1 00,300) GeV (M,,M;) = (1 00,500) GeV 




FIG. 3: Exclusions due to the AT?;, constraint, x%^, for 
(Ml, M2) = (100, 300) GeV, ^ = 0. Also shown is the region 
excluded by the B-physics constraints, and the 0% contour 
from Fig. [T] 



|i=200 GeV '^b 

(M,,M;) = (1 00.500) GeV (M,,M;) = (1 00,500) GeV 




1 1° tan/? 1 1° ioniS 

FIG. 4: ARb constraint, Xr^, for (Mi, M2) = (100, 300) GeV, 
fj, — 200 GeV. Also shown is the region excluded by the B- 
physics constraints and the 0% contour from the unitarity 
constraints. 



less constrained by the unitarity constraints, but for the 
higher value of M2 (right panel) the Rb modification that 
is caused by neutral-Higgs exchange starts to exclude 
high values of tan (3. 



B. LEP2 non-discovery 

The non-discovery of a neutral Higgs boson at LEP2 is 
relevant only for AIi < 114.4 GeV [62] ■ However, it is well 
known that these searches do not exclude certain other 
light neutral Higgs bosons, if they couple more weakly 
to the Z boson, or if they decay to final states that are 
more difficult to detect and identify. 

The constraint is implemented in an approximate way 
as follows. Following [6^ (see also [l^l), we consider the 
searches for a neutral Higgs boson that decays to bb jets or 
a tau pair. Thus, the quantity of interest is the product of 
the production cross section (proportional to the square 
of the ZZHi coupling) and the bb (or r+r") branching 



ratio (proportional to the square of the Hibb coupling). 
We denote the reduced sensitivity, as compared to the 
SM sensitivity, a "dilution factor" 63]: 

fT2HDM {ZHi -> Zbb) = asM {Zh Zbb) X {ZHi ^ Zbb) , 

(4.4) 

where the dilution factor is the product of a fac- 
tor C^{ZHi) related to the production and another, 
C^{Hi bb), related to the branching ratio: 

C^iZHi Zbb) = C^iZHi) X C^{Hi bb). (4.5) 

For the general 2llByLC'^ {Z Hi Zbb) = C|hdm is 
given by Eq. (4.3) in [llj. It depends on the neutral- 
sector rotation angles, as well as on tan/3, and is the 
same for bb and r+r^ final states. 



Ml 


80 GeV 


100 GeV 


114.4 GeV 


bb 


0.06 


0.25 


1.0 


r+r- 


0.06 


0.2 


1.0 



TABLE IV: Experimental suppression factors Cfxp [63 |. 



As experimental constraints, we approximate the 95% 
C.L. bounds obtained in (63| by linear interpolations 
passing through the points given in Table llVi and form 
the ad hoc single-sided penalty (summed over bb and 
T+r-): 

/ (^2 _(^2 \ 2 

2 c nnl ^2HDM ^oxp \ r.\ 

XLEP2 = 5.99 1^ — Y^Tcp ) (^-^^ 

for C|hdm > C'cxp and Mi < 114.4 GeV, and zero other- 
wise. The coefficient 5.99 corresponds to 95% probability 
in the context of our two "degrees of freedom" , tan /3 and 
Mh±. 



(M,,M;) = (1 00,300) GeV (M,,M;) = (1 00,500) GeV 




FIG. 5: LEP2 constraint, Xlep2- Left: (Mi,M2) = 
(100,300) GeV; Right: {Mi, Ah) = (100,500) GeV, n = 
in both cases. Also shown is the region excluded by the B- 
physics constraints, and the 0% contour from Fig. [1] 

We show in Fig. O the resulting exclusion for Mi = 
100 GeV, ^ = and two values of M2 as indicated. For 
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the higher value of M2, we note that some part of the 
otherwise allowed parameter space gets excluded. Given 
that C|jjj3M is determined by tan/3 and the a parame- 
ters, on might wonder why this excluded region depends 
also on Mh± ■ The reason is of course that the subspace 
of a that is allowed by the positivity and unitarity con- 
straints depends on Af^± , and need not overlap with the 
corresponding subspace of a for which CjjjDjyi is within 
the allowed range. 



> 

ID 

O 600 



400 



200 



/i=200 GeV 
(M,,M;) = (1 00,300) GeV 



LEP2 



ALLOWED 




400 



(M,,M;) = (1 00,500) GeV 




tan/3 




ALLOWED 
B physics 




ianfS 



FIG. 6: LEP2 constraint, Xlep2- Left: (Mi,M2) = 
(100,300) GeV; Right; (Mi,M2) = (100,500) GeV, fi = 
200 GeV in both cases. Also shown is the region excluded 
by the _B-physics constraints, and the 0% contour from the 
unitarity constraints. 



For = 200 GeV, the unitarity constraints no longer 
exclude high values of tan/3, but we see that the LEP 
non-discovery excludes high values of charged-Higgs mass 
(via its impact on the rotation matrix). In fact, at large 
values of tan/3, if Mi < 114.4 GeV, the LEP2 non- 
discovery requires (see Eq. (4.3) of [T]| ) 



1 



cos^ /3 



(4.7) 



For small values of cos/3, this is satisfied in three separate 
regions [see Eq. (|1.2p and Fig. [15] : 



(i) : ai ~ 0, 

(n) : \a2\ — 7r/2, 

(iii) : |ai|~7r/2 and a2 

C. The p parameter 

The p parameter, defined as 



(4.8) 



Ml cos2 6'w ' 



(4.9) 



is very sensitive to fields that couple to the W and Z [65 
Experimentally, p is constrained as [ssj 



p^^p = 1.0050 ±0.0010. 



(4.10) 



The deviation from unity is mostly due to the top-quark 
one-loop contributions, but there is also a weak depen- 
dence on the SM Higgs mass. In order to extract this 
quantity from the data, one fits for a SM Higgs mass, the 
contribution of which should then be subtracted from the 
2IIDM prediction. 

The additional Higgs fields of the 2HDM can easily 
spoil the agreement with the SM 7]. Roughly speaking, 
this constraint requires the Higgs masses to be not too far 
from the W and Z masses, and not very much apart. For 
the general CP-non-conserving 2HDM, the results of 
for the contributions to p, denoted Ap — /92HDM — Psm, 
were generalized in [llj and given by Eqs. (4.8)-(4.12) 
there. 



> 

O 600 



400 



200 



/x=0 

(M,,M;) = (1 00,500) GeV 



Ap 




200 



(M,,M;) = (1 00,500) GeV 



tan^ 




tan(5 



FIG. 7; Exclusions due to the Ap constraint, Xpi for 
(Mi,M2) = (100,300) GeV, and (A/i,M2) = (100,500) GeV, 
p = in both cases. Also shown is the region excluded by the 
-B-physics constraints, and the 0% contour from Fig. [T] 



In order to study this constraint, we evaluate 



defined by 



'2HDM 



(4.11) 



with a{p) — 0.0010, subtracting the SM contribution cor- 
responding to an SM Higgs mass Mq = 129 GeV. 

In Figs. [7] and [8] we show the impact of xj, in constrain- 
ing the parameter space for /x = and p = 200 GeV, 
respectively. At the higher value of M2, and high values 
of tan /3, this constraint tends to exclude both low and 
high values of M^± . The irregular boundaries seen in 
Fig. [8] are obviously due to the scanning not finding the 
whole allowed region. (At high p and high tan/3, only 
very tiny regions in ai and a2 are allowed [19].) 

It is instructive to consider the contribution to Ap 
in the limit of large values of tan/3. The expressions 
(4.10) and (4.12) of 11] simplify considerably in the limit 
sin /3 — > 1 , and provide an understanding of features seen 
in Figs. [7] and 11 

<w(O)-cos2 0H'Aff(O) 



647r2 ^ 
j 
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> 

O 600 



400 



200 



/i=200 GeV 
(M,,M;) = (1 00,300) GeV 



Ap 




600 



400 



(M,,M;) = (1 00,500) GeV 



tan/3 




tan(5 



FIG. 8: Exclusions due to the Ap constraint, Xpi for 
(Ml, Ma) = (100,300) GeV, and (Mi,M2) = (100,500) GeV, 
/i = 200 GeV in both cases. Also shown is the region excluded 
by the B-physics constraints, and the 0% contour from the 
unitarity constraint. 



J2iRjiRk3 - RkiRjsfFApiMf, M^)] , (4.12) 



and 



<'^(O)-cos2 0v^Aff(O) 



Gin 



: [E (3^Ap(M|, Mf) - SFaMw, Mf) 

j 

+ 3i^Ap(M^, M^) - 3FAp{Ml M^)] (4.13) 

with FAp{ml,ml) found in Eq. (4.11) of ll]. This func- 
tion vanishes when the two masses are equal, and grows 
quadratically with the bigger of the two masses. The 
contribution given by (|4.13[) is rather small, since Mz 
and Mw are relatively close. 

The only part sensitive to the charged-Higgs mass is 
the first sum in (|4.12p , which may however get significant 
contributions from all three neutral Higgs bosons, j = 
1,2,3. Consider first the contribution from Hi. Barring 
cancellations (see below), a viable parameter point must 
for large tan (3 and large ^ Mi have 



Rli+Rl,^l. 



(4.14) 



Expressed in terms of the angles, this means that cfc^ 

2 
2 



S2 < 1, or 



7r/2, and a2 ~ 0. 



(4.15) 



Note that this condition is compatible with (|4.7p and 
corresponds to case (in) of (|4.8p . 

Next, we consider the case (with increasing splitting 
between M2 and M3, the contribution to Ap will in- 
crease) 



then again barring cancellations, the Ap condition re- 
quires 



Y^iRh + R%) = + sl + cl « 1. 



(4.17) 



j=2 



This condition cannot be satisfied, and low values of M[j± 
are thus excluded when tan /3 is large and M2 is large. 

For high values of M[j± , the most sensitive contribu- 
tion is that involving Hi. But unless Mh± is close to 
M2, the condition (|4.15p . related to the Hi contribution, 
is not sufficient. The H2 and H^ contributions are pro- 
portional to i?2i + ^23 ^-iid R\i -\- i?§3, respectively. In 
the limit (|4.15p . these are both equal to 1, i.e., when 
the rotation matrix is adjusted such as to cancel the Hi 
contribution to the first sum in (j4.12p . there is no sup- 
pression of the H2 and H-^ contributions. As a result, 
such high values of M[j± are forbidden. 

An exception to this situation arises when p, is large 
compared with M2. Then, there can be a considerable 
splitting between M2 and M3, Mu± and Afa can be sim- 
ilar, and a cancellation between the (Mi, Mu±) and the 
(Ml, Afa) terms of (|4.12p is possible. As a result, for 
large values of /i, large values of Mh± can be allowed. 

Finally, we consider the situation 



Afi < M2 ~ Afa ~ Mh± , 



(4.18) 



Mh± < Af2 ^ M3 



(4.16) 



which applies to the "finger" protruding to the right in 
the right panel of Fig. [H Then, the contributions of the 
two sums in Eq. (|4.12p tend to cancel. For j = 1 and 

fc = 2, 3, we get 

[i?ii + i?i3 — (i?iii?23 — -^21^13)^ 
- (i?lli?33 - R3lRl3f]FApiMl Mi) = 0, (4.19) 

where we have used the orthogonality of the rotation ma- 
trix. Thus, when tan/3 is large, and Af2 — A^3 ~ A^^±, 
there is a cancellation among the terms in (|4.12p . only 
the (small) (|4.13p part of the Ap constraint is relevant. 



D. Muon anomalous magnetic moment a^j 

The precisely measured muon anomalous magnetic mo- 
ment |66| . 

a^-exp = 5(5-2)p,exp = 11659208(5.4)(3.3) x 10-^" 

(4.20) 

is a sensitive probe of new physics (for a recent review, 
see [sll). The statistical and systematic uncertainties 
(given in parentheses) combine to an over-all uncertainty 
of 6.3 X 10^^". The corresponding SM prediction, includ- 
ing weak and strong effects, is [69| 

ap,sM = 11659179.3(6.8) x 10"^°, (4.21) 

creating a 3f7 "tension" with the experimental result. 
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In the 2HDM, there are additional contributions, dom- 
inated by the two-loop Barr-Zee effect (T^I with a pho- 
ton and a Higgs field connected to a heavy- fermion loop. 
For the CP-conserving case, the contribution is given 
by [i [Zii- For the general (CP-violating) 2HDM, the 
top-quark contributions to for the muon, is given by 
Eq. (4.13) in [ll|, whereas the 6-quark contributions to 
the fermion loop is given by 



tan^ /3i?,3 



1 



R2 f 

cos2/3 ^■1-' ^AfJ 



(4.22) 



with A^c = 3 the number of colours associated with the 
fermion loop, ac.m. the electromagnetic finestructure con- 
stant, Qh = —1/3 and mf, the 6-quark charge and mass, 
and the muon mass. The r-loop contribution, which 
we also include, is given by a similar expression, with 
obvious substitutions for the colour factor, charge and 
mass. The functions / and g are given in [TOj . 

At low values of tan/3, these contributions are negli- 
gible, but the h- and r-loop contributions can become 
relevant at very large values of tan (3. As a measure of 
the possible conflict with the 2HDM, we consider 



Aa,, 



Aa 



;i,2HDM 

cr(ap) 



(4.23) 



where Aa^_2HDM is the 2HDM-specific contribution, and 
for the uncertainty we take the (SM) theoretical value, 
(j{a^) = 6.8 X 10^^*^, since this is larger than the experi- 
mental one. 

Actually, this constraint does not have any significant 
impact within the range of tan /3 considered. Let us con- 
sider its "natural value" as that contributed by one of the 
two terms in (|4.22p . with the rotation matrix element set 
to 1. This reaches = ^^(1) for tan/3 of the order of 
70. However, there can be significant reductions by the 
rotation matrix elements, and also cancellations among 
the two terms. At such high values of tan/3 the B — > tVt 
constraint (see Sec. IIII B[) is important at low values of 
Mjj± , and the unitarity constraint may be important at 
high Mh± (depending on the relative magnitude of M2 
and /i). 

We recall that at high tan/3, the rotation matrix is 
rather constrained by unitarity [l9| . Let us focus on 
the contribution of the lightest neutral Higgs boson, 
whose contributions to (|4.22p are given by R\i and R\^. 
Thus, in spite of the enhancement given by the tan/3- 
dependent factors, this contribution to may be small 
if 



\ai \ ^ 7r/2, and a-i ^ 0. 



(4.24) 



In the limit of large tan /3 and large Mh± , this is actu- 
ally the only region allowed by the Ap constraint, see 
Eq. (|4.15p . We conclude that at large tan (3 and large 



Mh± , the constraint is covered by the Ap constraint. 
But this coefficient Rf^ + R\^ arises in one case from 
the Yukawa couplings, and in the other from the gauge- 
Higgs couplings. Furthermore, at large tan/? and mod- 
erate Mu± values, the constraint is covered by the 
B TVr constraint. 



E. Summary on neutral-sector constraints 

We will here summarize the conclusions on the neutral- 
sector constraints, treating first the simpler case of /j, = 
(where tan /? is bounded) , and next comment on the less 
restrictive case of "large" p (where also larger values of 
tan /3 are allowed) . The Rh constraint is at low values of 
tan /3 dominated by the charged-Higgs-exchange contri- 
bution. This part of the Rh constraint is thus indepen- 
dent of the neutral sector. 



1. The case /i = 

For M = 0, and (Mi,Af2) = (100,300) GeV, the only 
neutral-sector constraint that has some impact, apart 
from Rh at low values of tan/3, is the Ap constraint, 
which excludes the higher range of Mh± , as illustrated 
in the left panel of Fig. [T] However, for M2 = 500 GeV, 
other constraints are also important. The LEP2 non- 
discovery rules out large values of tan /3 and Mfj± , and 
to some extent, also the Ap constraint rules out some 
region of large tan/3. 

However, with p = 0, high values of tan [3 are also ex- 
cluded by the unitarity constraints, and, to some extent, 
the low values of Mfj± are excluded by the B-physics 
constraints. 



2. The case p > Mi 

When p > Ml , large values of tan (3 become accessible. 
This parameter region is known as the decoupling region 
72]. We here distinguish two cases 



{i) Ml < p < M2, 
{a) Ml < M2 < p. 



(4.25) 



The Rh constraint, which at low tan/3 is dominated by 
the charged-Higgs contributions, can at large tan/? also 
exclude some region of neutral Higgs boson mass values 
(compare the left and right panels of Fig. 0]). Further- 
more, the LEP2 constraint may exclude high values of 
Mfj± (see right panel of Fig. O and the Ap constraint 
may at high tan (3 constrain the range of M^± values to 
a band around M2. 
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V. COMBINING ALL CONSTRAINTS 

Let us now combine all constraints. This is done by a 
dedicated scan over d for each point in tan (3 and M^± . 
The value of is determined as 



\ Xgcncral ' 



(5.1) 



where Xgenorai is given by Eq. (|3.19|) and the sum runs 
over the observables LEP2 non-discovery, Ap and 
a^, all of them evaluated at the same point in a. This 
quantity is then minimized over a: 



= min x^, 



(5.2) 



for fixed tan [3 and Mu± and allowed regions are deter- 
mined. In the tan/3-M^± plane, the allowed regions will 
in general be less that the intersection of the regions that 
are allowed by the individual constraints. The reason 
is that the individual constraints may refer to different 
parts of the three-dimensional a space. 

We shall split this discussion into the two cases /i = 
and Ml < fj,. In the former case, unitarity restricts the 
allowed range of tan/3, as illustrated in Fig. [l] whereas 
in the latter case also higher values of tan /3 are allowed. 



A. Combining all constraints for ^ = 

In Fig. M we display the 90 and 95% CL. hmits 
for (Mi,M2) = (100,300) GeV and (100,500) GeV, in 
both cases for — 0. For the case of moderately low 
M2 = 300 GeV, we note that there is little additional ex- 
clusion, other than that due to the i?-physics constraints 
and unitarity. The little extra is due to the Ap constraint, 
at high values of Af^± . However, for A/2 = 500 GeV, 
there is a considerable reduction of the allowed parame- 
ter space at tan > 1 — 1.5. In this range of tan f3 values, 
we see from Fig. [S] (right panel) , that high values of M^± 
are excluded by the LEP2 non-discovery constraint, and 
from Fig. [7] (right panel), we see that low values of Af^± 
are excluded by the Ap constraint. 

Similarly, Fig.[in]is devoted to the case Mi = 150 GeV. 
In this case, the LEP2 non-discovery plays no role, 
but there are of course also other differences, due to 
the way different parameters are correlated by the con- 
straints. Overall, this case is less constrained than the 
Ml = 100 GeV case of Fig. M 



B. Combining all constraints for A/i < ^ 

When Ml < /i, the unitarity constraints no longer re- 
strict tan P to low and moderate values. However, we 
shall see that various other constraints may cause a cut- 
off for large tan/3. For Mi = 100 GeV and two val- 
ues of Af2, namely 300 GeV and 500 GeV, we display in 



/j.=0 

(M,,M;) = (1 00,300) GeV 



ALL 




(M,,M^) = (1 00,500) GeV 



600 



400 



tan|8 




tan(5 



FIG. 9: Exclusions due to all constraints, for (Afi,M2) = 
(100,300) GeV, and (Mi,M2) = (100,500) GeV, both with 
fj. = 0. Heavy dashed: 90% C.L., heavy solid: 95% C.L. Also 
shown is the region excluded by the B-physics constraints 
(thin dashed and solid), and the 0% contour from Fig. [T] (thin 
solid) . 
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(M,,M;) = (1 50,500) GeV 



10 



tan|8 




tanjS 



FIG. 10: Similar to Fig. M for Mi = 150 GeV. 



Figs. fTT] - [T3l the allowed regions for a range of p-values, 
from 200 GeV to 600 GeV. 

For the lower value, M2 = 300 GeV (left panels), the 
allowed region in the tan /3-Mjy± plane is fairly extended, 
whereas for the higher value, M2 = 500 GeV, it is more 
constrained, until fi reaches values comparable to M2 (see 
the right panels of Figs. [IH [H] and [H]) . 

For /i = 500 GeV and M2 = 500 GeV (not shown), the 
exclusion of low and high values of Mh± for tan/3 > 2 
is due to the Ap constraint, whereas the LEP2 non- 
discovery by itself does not exclude anything in this case. 
However, the simultaneous imposition of the LEP2 non- 
discovery and the Ap constraints yields a "forbidden fin- 
ger" at tan/3 ~ 1.5 and low Mij^± , as well as some exclu- 
sion at low tan /3 and high Mff± . 

In all cases, we note that simultaneously high val- 
ues of both tan/3 and Mh± are excluded, except when 
M2 < /i. This is due to the Ap constraint, as discussed in 
Sec. HTCl For the cases shown (Mi = 100 GeV), also the 
LEP2 non-discovery plays a role, but this will of course 
not have any impact for Mi > 114.4 GeV. 
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/i=200 GeV 
(M,,M;) = (1 00,300) GeV 



ALL 




600 



400 



(M,,M;) = (1 00,500) GeV 




tan/3 



FIG. 11: Similar to Fig.[9l for fi = 200. 



tan(5 



A. Low values of tan f3 

At low values of tan/?, both low and high values of 
fi lead to consistent solutions, with values of ai and a2 
distributed over extended regions of the parameter space, 
as shown in Fig. [M] for the case fi = 400 GeV, Mi = 
100 GeV and two values of M2, namely 300 GeV and 
500 GeV. (Further plots of this kind, but taking into 
account only the positivity and unitarity constraints, are 
presented in [Toj.) 

We have here plotted the distributions of all a\ and 
a2 for which the total ^ < 5.99 (see Eq. (|5.ip ). i.e., 
in general several points in (q;i,q;2) for each point in 
{t&n (3, Mh± )- 



> 

O 600 



400 



200 



/i=400 GeV 

(M,,M;) = (1 00,300) GeV 



ALL 




500 



400 



(M,,M;) = (1 00,500) GeV 




tan^ 



1 1° tan/3 

FIG. 12: Similar to Fig. [H for = 400 



VI. PROFILE OF SURVIVING PARAMETER 
SPACE 



7t/2 



tan/3 < 5 
M,,M2,^) = (100,300,400) 



(100,500,400) 



-7t/2 





i 








If 


If ii 


— 







-7t/2 



7t/2 -7t/2 



7t/2 
a, 



FIG. 14: Normalized distributions of allowed regions in the 
a\-cx2 space, for low values of tan/3 values. Contours are 
shown at each negative power of 10, as appropriate. Yel- 
low (light blue) indicates where the normalized distribution 
is higher than 10"* (3 x 10~^); green (purple) levels above 
10"^ (3 X 10"^). Along the lines 02 = ±7r/2, Rx is CP-odd, 
and there is no CP violation. 



We shall here give a profile of the surviving parameter 
space in terms of three "hidden" parameters, ai, and 
M3. In discussing the surviving parameter space, we shall 
distinguish between low and high values of tan /3, giving 
some details relevant to tan/3 < 5 and tan/J > 10. 



/i=600 GeV 

(M,,M;) = (1 00,300) GeV (M,,M^) = (1 00,500) GeV 




FIG. 13: Similar to Fig. [H for = 600. 



At very low tan /?, an important constraint is the pos- 
itivity of A2: 

A2 = -^{s\c\mI -f (ciC3 - SxSis-ifMl 

+ (C1S3 + s^S2CzfMl - cli?\ > 0, (6.1) 

and the constraint from unitarity that it does not become 
"large". For small /i, the M|-term cannot be too large 
(in order not to violate unitarity). This requires 

IC1S3 + S1S2C3I < 1. (6.2) 

As fi becomes large, with M2 fixed, the M|-term must 
compensate the /i^-term, with |ciS3 -I- S1S2C3I = 0{1). 
The distributions in Fig. [14] are seen to satisfy this con- 
dition. 

While the allowed range of Mh± depends on the neu- 
tral Higgs boson mass M2, it is typically of the order of 
300-700 GeV. In most of the allowed parameter space, 
CP is violated, but along the edges a2 —> ±7r/2 there 
is no CP-violation. This is the limit where the lightest 
Higgs boson, Hi, becomes CP odd. 
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B. High values of tan/3 

We show in Fig. [15] the populated regions in the ai- 
a2 plane, for /i — 400 GeV and tan/? > 10. Two points 
are worth noting: (1) The allowed regions satisfy the con- 
straints of (|4.8p . (2) The majority of points do not satisfy 
the condition (|4.15p . meaning that the Af^± ^ Mi case 
is not very relevant here. Instead, the degenerate case, 
/i ~ Mh± ^ M2 ^ is important, as also reflected in 
Table El c). 



tan/? >10 

(M,,M2,^) = (1 00,300,400) (1 00,500,400) 
7t/2| 



-7t/2 








a, 



7t/2 — 7T/ 



7t/2 



FIG. 15: Normalized distributions of allowed regions in the 
ai-a2 space, for high tan /3 values. Contours and colour codes 
are as in Fig. 1141 with additionally red (blue) above 10~^ 
(3 X 10"^). 



be low, just marginally above M2, as seen in Table IVl a) . 
This pattern is valid also for Mi < /j, provided only that 
/i < M2, see Table fVl b) and c). For M2 < /i, on the 
other hand, M3 can be large. 

a) = 0, (Afi, M2) = (100, 300 [500]) GeV 



tan f3 


e< 1.1 


l.K ^ < 1.5 


1.5 < C 


5-10 
< 5 


94.6 [ 0.0]% 
43.4 [88.4]% 


5.2 [ 0.0]% 
49.8 [11.6]% 


0.2 [ 0.0]% 
6.8 [ 0.0]% 


b) /i = 200 GeV, (Mi, M2) = (100, 300 [500]) GeV 


tan/3 


^< 1.1 


l.K C < 1-5 


1.5 < C 


> 10 

5-10 
< 5 


74.0 [95.8]% 
49.0 [91.4]% 
30.9 [81.5]% 


24.9 [ 4.2]% 
48.4 [ 8.6]% 
56.8 [18.5]% 


1.0 [ 0.0]% 
2.6 [ 0.0]% 
12.3 [ 0.0]% 


c) /i = 400 GeV, (Ml, M2) = (100, 300 [500]) GeV 


tan /3 


e< 1.1 


l.K C < 1-5 


1.5 < C 


> 10 
5-10 
< 5 


0.0 [91.6]% 
0.0 [76.3]% 
0.0 [64.1]% 


69.0 [ 8.4]% 
47.8 [23.7]% 
15.3 [35.8]% 


31.0 [ 0.0]% 
52.2 [ 0.0]% 
84.7 [ 0.2]% 


d) /i = 600 GeV, (Afi, M2) = (100, 300 [500]) GeV 


tan/3 


C< 1.1 


l.K e < 1-5 


1.5 < ^ 


> 10 
5-10 
< 5 


0.0 [ 0.0]% 
0.0 [ 0.0]% 
0.0 [ 0.0]% 


0.0 [93.8]% 
0.0 [92.6]% 
0.0 [94.5]% 


100.0 [ 6.2]% 
100.0 [ 7.5]% 
100.0 [ 5.6]% 



TABLE V: Distribution of M3 values, £, = M3/M2 



At high values of tan /3, /i has to be comparable to 
M2, or higher. In particular, no solution exists for high 
values of tan/3 and ^ = 0. This is due to the unitarity 
constraints. Also, we note from Figs. [9HT31 that unless /i 
is comparable to, or larger than M2, the allowed range in 
tan (3 and Mh± can be rather limited. The range of Mh± 
depends on the neutral Higgs boson mass AIi , being typi- 
cally of the order of 400-600 GeV. In order to understand 
the large-tan (3 parameter space, let us review 



Ai 



1 



[clcjM^ + (C1S2S3 + sic^fM^ 

2 ?\/r2 



+ (C1S2C3 - sis^ifMi 



s}fi''] > 0. 



(6.3) 



This should be positive, but not "too large" . High values 
of /i require high M3 and |ciS2C3 — S1S3I = 



C. Distribution of M3 

We recall that with our choice of parameters, the third 
neutral Higgs mass is a derived quantity. The distribu- 
tion of M3 can be discussed in terms of the dimcnsionless 
ratio 



M3 

M2' 



i<e 



(6.4) 



For three bins in ^, this is for /i = 0, 200, 400 and 600 GeV 
distributed as given in Table IVl For fj, = 0, M3 tends to 



D. The Standard-Model-Uke limit 

The parameters of the 2HDM can be chosen such that 
the ZZHi, bbHi and ttHi couplings all approach the 
corresponding SM values. This requires, in our notation 



cos(/3 — ai) cos 0:2 
cos ai cos a2 

cos/3 
sinai cos 0:2 
sin/3 



1, 
1, 

1, 



(6.5) 



which is satisfied for /3 ~ ai and 0:2 — 0. There could 
also be a "quasi-SM-like" limit, where one or more of the 
above quantities approaches —1 (denoted "Solution B" 
in [zl]). For the familiar observables, such a sign change 
would not have any effect. 

For tan/3 ^ 5, the allowed regions in the 01-0:2 space 
are rather extended, and SM-like solutions are found for 
a range of mass values. For tan (3 > 10, on the other 
hand, the populated parts of the 0:1-02 space become 
very localized, and have the following features: (i) for 
M2 ^ /i, oi ~ and [02 1 > 0, and (ii) for M2 ^ /i, 
additional regions emerge for small values of ]o2| and 
]oi] ~ 7r/2. The latter, seen as small specks near the 
horizontal axis in the left panel of Fig. [151 correspond to 
the SM-like (and "quasi-SM-hke") case. 



15 



For the case Mi = 100 GeV, which is studied in most 
of our figures, there is of course a Umit to how close 
we come to the SM hmit, since an SM Higgs mass of 
this low value is excluded. Actually, the fact that for 
(Ml, M2) = (100, 500) GeV and fJ. = 0, some region of the 
tan/3-MjLf± plane is excluded by the LEP2 non-discovery 
constraint (see right panel of Fig. [6|), means that those 
regions correspond to solutions near the SM limit. 



VII. POSSIBLE FUTURE CONSTRAINTS 

It is interesting to see how the parameter-space con- 
straints would be modified by possible future results from 
the B-physics sector. We shall not consider any change 
to the Ap or LEP2 constraints, but rather discuss the 
possibility that the central value for the branching ratio 
for B be reduced from 3.55 x 10"'' to 3.20 x 10"'', 

a value closer to the SM prediction. Also changing the 
overall uncertainty to 0.25 x 10^^, we find that the result- 
ing constraints are significantly modified, as illustrated in 
Figs. [IS and [17] for = and 400 GeV, respectively. 
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FIG. 16: Similar to Fig. [H for BR{B 
and an over-all uncertainty a[B{B —* 
(c.f. Eq. ((33J). 



^.7)] = 



tan^ 
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FIG. 17: Similar to Fig. [H for /i = 400 GeV. 



VIII. SUMMARY 



We have shown that the B-physics results, together 
with the precise measurement of the p-parameter at LEP 
and the constraint of tree-level unitarity of Higgs-Higgs 
scattering, exclude large regions of the 2HDM (II) pa- 
rameter space. High values of tan (3 are excluded unless 
both M2 and M3 are heavy. Furthermore, they should 
be reasonably close to each other. Improved precision 
of the B — > measurement could significantly reduce 
the remaining part of the parameter space, but it appears 
unlikely that the model could be excluded other than by 
a negative search at the LHC. 

What is the corresponding situation for supersymmet- 
ric models? While the consistency then is guaranteed by 
internal relations, it should be kept in mind that light 
charged Higgs bosons would be in conflict with the B- 
physics data unless some superpartner (for example, the 
chargino (tI]) is also light. A possibility which has re- 
ceived some attention, is a light chargino and a light 
stop [zHl. It has also been shown that anomalous ef- 
fects at large tan/3 could weaken the bound on Mh± 
without light superpartners [t^ (see also [z3)- However, 
the more recent data on _B ^ tVt discussed in Sec. IIII BI 
would presumably close this loophole (see Fig. 2 of 17^). 

Similar scans over the parameter space of the Con- 
strained MSSM [7i,[83,[8ll differ from the present work in 
one major respect: they are required to yield an amount 
of dark matter that is compatible with the WMAP data 
[s^ . Additionally, the constraint is more severe, due 
to one-loop contributions involving superpartners of the 
muon and muon neutrino. These studies are also more 
focused on the high-scale parameters, like mo and mi/2, 
with less emphasis on M^± and tan (3. The study by 
[sol shows a preference for positive ^ (higgsino mass pa- 
rameter, not to be confused with the /j, of Eq. (|2.ip ). a 
relatively light charged Higgs mass (a few hundred GeV) 
and rather high values of tan /3 (~ 50 — 60) . 
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